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TRANSPORT  VIA  MOMENTS  OF  QUANTUM  DISTRIBUTION  FUNCTIONS 
H.  L.  GRUBIN  and  T.  R.  GOVINDAN 


INTRODUCTION 


The  development  of  crystal  growth  techniques  has  made  it  possible  to  fabricate  devices  with 
very  small  dimensions  and  sharp  interfaces.  Consequently,  a  variety  of  device  concepts 
based  upon  small  dimensions  with  sharp  interfaces  has  emerged.  These  include  such  high 
speed  devices  as  the  HEMT,  HBT,  quantum  resonant  tunneling  devices,  quantum  wire 
structures  and  quantum  dots.  The  structure  and  operation  of  these  devices  is  based  upon  a 
perceived  physical  picture,  often  based  upon  results  of  quantum  transport  theory. 


While  all  devices  are  governed  and  sometimes  limited  by  quantum  mechanical  properties, 
many  devices  do  not  require  quantum  mechanics  for  an  understanding  of  their  basic 
operation,  but  do  need  quantum  mechanics  to  provide  key  electrical  properties.  For 
example,  heterostructure  bipolar  transistors  sustain  low  levels  of  current  at  low  bias  levels; 
these  currents  are  dominantly  tunneling  currents.  Thermionic  contributions  to  current 
occur  at  high  bias  levels.  Until  recently,  the  drift  and  diffusion  equations  as  well  as  the 
moments  of  the  Boltzmann  transport  equation  did  not  include  a  description  of  tunneling 
currents  in  any  way  other  than  through  ad  hoc  arguments.  Tunneling  requires  quantum 
mechanical  contributions. 


The  recent  advance  in  the  description  of  transport  in  devices  with  quantum  contributions  is 
due  to  the  discussion  of  Ancona  and  Iafrate  (1989),  who  demonstrated  that  the  quantum 
corrections  to  multi-particle  transport,  first  discussed  by  Wigner  (1932),  could  be 
incorporated  in  the  drift  and  diffusion  equations  by  adding  to  the  classical  potential  a  term 
proportional  to  a  quantum  mechanical  based  potential: 

(1) 

where  the  subscripts  denote  derivatives.  Within  the  framework  of  the  drift  and  diffusion 
equation  the  ’Q’  term  modifies  the  particle  current  as  follows: 

(2)  j(x,t)=n^kbT[(V  +  aQ)/kbT  +  in(n)]x 

The  proportionality  constant  ’ a’  was  evaluated  [Ancona  and  Iafrate  (1989)]  in  the  high 
temperature  limit  from  the  equation  of  motion  of  the  Wigner  function  and  was  found  to  be 
equal  to  1/3.  The  factor  T/3’  is  a  serious  issue,  in  that  intuitively,  as  well  as  through 
arguments  associated  with  the  single  particle  Schrodinger  equation,  it  would  be  anticipated 
that  the  constant  V  would  be  unity. 

The  problem  addressed  in  the  last  paragraph,  namely  the  origin  of  the  factor  a = 1/3,  as  well 
as  the  procedures  under  which  quantum  contributions  to  the  balance  equations  are 
obtained  formed  the  basis  of  the  study  under  Contract  DAAL03-90-C-0005.  This  was 
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complemented  by  numerical  simulations,  using  ’characteristic’  algorithms  to  examine  the 
extent  to  which  the  quantum  corrected  results  are  relevant. 

The  means  by  which  the  problem  is  addressed  is  through  an  expansion  of  the  equation  of 
motion  of  the  Wigner  function,  and  an  expansion  of  the  equation  of  motion  of  the  density 
matrix.  The  discussion  of  the  Wigner  function  is  in  the  form  of  a  review,  in  that  a 
comprehensive  treatment  appears  in  Grubin  and  Kreskovsky  (1989).  The  work  on  the 
density  matrix  is  new. 

QUANTUM  CORRECTED  EQUATIONS  OF  MOTION  FOR 
THE  WIGNER  FUNCTION  AND  THE  DENSITY  MATRIX 

The  Wigner  Equation  of  Motion:  The  analysis  under  this  study  was  restricted  to  one  space 
and  one  momentum  direction  and  includes  Fokker-Planck  scattering.  The  equation  of 
motion  of  the  Wigner  function  with  Fokker-Planck  dissipation  as  discussed  by  Stroscio 
(1986)  is: 

(3)  ft  +  (p/m)fx + (1/ifi  )(l/2*ft)_J  +  -  dp’^/  +  -  dx’f(p’pc)[V(xpc’)t)]exp[i(p-p’)x,M  ] 

=27[pf]p+Dfpp 


where: 

(4)  V(xpc\t)  =  [V(x +xV2,t)-V(x-xV2,t)] 

The  potential  energy  V(x,x’,t)  is  referred  to  below  as  the  Liouville  potential.  As  discussed 
by  Frensley  (1990)  the  first  term  of  the  Fokker-Planck  dissipation  corresponds  to  a 
frictional  damping  term,  the  second  corresponds  to  thermal  fluctuations.  This  will  be 
discussed  in  more  detail  below. 

It  is  direct,  but  nontrivial  to  demonstrate  that  the  integral  in  equation  (3)  reduces  in  the 
classical  case  to  Vxfp*,  thus  equation  (3)  reduces  to  the  Boltzmann  transport  equation  with 
Fokker-Planck  dissipation: 

(5)  ft  +  (p/m)fx-Vxfp=2-r[pf]p  +  Dfpp 

To  second  order  in  n ,  the  Wigner  equation  of  motion,  or  equivalently,  the  quantum 
corrected  Boltzmann  equation  is: 

(6)  ft + (p/m)fx-Vxfp  +  {h2  /24)Vxxxfppp = 27  [pf]p + Dfpp 

The  left  hand  side  of  equation  (6)  has  been  discussed  in  detail  by  Ancona  and  Iafrate 
(1989)  and  Grubin  and  Kreskovsky  (1989). 

The  Equation  of  Motion  of  the  Density  Matrix:  This  equation  has  been  discussed  in  the 
context  of  quantum  structures  by  Frensley  (1990).  While  the  equation  of  motion  is  a 
fundamental  equation,  it  is  equivalent  to  and  can  be  obtained  directly  from  the  Wigner 
equation  of  motion.  In  the  latter  case  multiply  equation  (3)  by  dp[exp[i2f  p/ft]]  and 
integrate  over  p.  For  this  integration  p  varies  from  -«>  to  +»,  ana  it  is  assumed  that  the 
Wigner  function  and  all  necessary  derivatives  with  respect  to  momentum  vanish  as  p=*±®. 
Identifying  the  density  matrix  through  the  Weyl  transformation: 

(7)  p  (x + f  ,x-r  ) = [l/(2itft  )]-o/0°dpfw(p^)exp[2ipf  In  ] 


the  transformed  equation  becomes  the  density  matrix  in  the  coordinate  representation: 

(8)  p  t  +  (ft /2mi)p x  -( 1/ift )  [  V(x  +  f ,t)- V (x  r,t)]p+27fpc  +  [4Df  2/ft  2]p  =  0 
An  expansion  of  the  Liouville  potential  about  the  point  x  yields: 

(9)  pt+(ft/2mi)pxr(l/ift)[2cVx  +  (f3/3)V^]p+27Cpf +[4Df2/ft2]p=0 
which  is  the  transformed  equivalent  of  equation  (6). 

Retaining  only  theJerms  that  are  linear  in  C  yields  the  transformed  equivalents  the 
classical  Boltzmann  equation.  As  will  be  seen  the  presence  of  the  factor  *1/3*  in  equation  (9) 
is  the  key  ingredient  in  obtaining  the  multiplicative  factor  in  the  quantum  potential.  It  is 
emphasized  that  quantum  corrections  to  0(ft 2 )  are  contained  in  the  term  [(C 3  /3) V^lp . 

Equations  (6)  and  (9)  are  the  relevant  equations  for  the  analytical  part  of  this  study. 
Equations  (3)  and  (8)  are,  however,  the  ones  to  deai  with  in  a  fundamental  approach  to 
quantum  transport.  Results  with  equation  (8)  are  introduced  later  into  the  discussion. 

THE  APPROXIMATE  EQUILIBRIUM  DISTRIBUTION  FUNCTION 

In  the  absence  of  dissipation  the  approximate  Wigner  distribution  function  to  second  order 
in  ft  is  Wigner  (1932): 

(10)  fw  =  exp-^[p2/2m4'U(x)]{l-(>  2^/4)[(Vxx-^Vx2/3)-^(p2/3m)Vxx]} 

Where  \  2  =  ft  20/2m,  and  p  =  l^k^T).  The  equivalent  approximate  density  matrix  solution 
to  order  0(ft  2 )  is: 

(11)  p(x-f-f  pc-f  )=  Nexp-[f 2  A  2  +£V]{1-(X  2^/6)(l  +  f  2A  2)Vxx  +  >  2/92  Vx2/12} 
where, 

(12)  N  =  2(m/2ir/?«2)1/2. 

Note:  First  the  classical  density  matrix  equivalent  to  the  Boltzmann  distribution  is: 

(13)  p  (x  +  f  ,x-f )  =  Nexp-[f  2  A 2  +  P V]; 

Second,  the  brackets  { •  •  • }  contains  quantum  corrections  to  the  density  matrix.  Third,  the 
approximate  density  matrix  is  real,  indicating  that  the  current  is  zero.  Fourth,  the  diagonal 
components  of  the  approximate  density  matrix,  which  are  obtained  for  f  =0,  are  equal  to: 

(14)  P(x*)  =  N[exp-/5V][l-(X  ^ 

This  result  also  emerges  from  an  integration  of  the  Wigner  function  over  all  values  of 
momentum  (Ancona  and  Iafrate  (1989). 

The  significance  of  equation  (14)  for  device  applications  is  that  it  describes  the  way  in 
which  quantum  corrections  alter  the  built-in  potential.  For  example,  when  it  is  recalled  that, 
classically,  the  potential  and  carrier  density  are  related  through  the  expression: 

(15)  p  =Nexp-0V 
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and  therefore  that  £Vx  =  -px/p;  then  through  incorporation  of  the  quantum  potential 
(equation  (1)),  it  is  straightforward  to  show  that  Q  =  (h  2/4m)^[Vxx-^Vx2/2],  and  that: 

(16)  p  =  N[exp-0V][l-(£Q/3)] 

Several  points  are  clear:  First  there  may  be  an  alteration  of  the  built-in  potential  arising 
from  gradients  in  the  carrier  density;  second,  the  factor  of  *3’  that  appears  in  equation  (16) 
is  seep  to  be  a  consequence  of  the  truncation  of  the  Taylor  expansion  of  the  potential  in 
both  the  density  matrix  and  Wigner  equations. 

Quantum  Corrected  Energy  Density  Matrix:  There  is  a  characteristic  energy  associated  with 
the  equilibrium  system  obtained  from  both  the  Wigner  function  and  the  density  matrix.  To 
obtain  it  from  the  density  matrix  we  define  an  energy  density  matrix  through  a  Weyl-like 
transformation: 

(17)  E(x  +  C  tX-C )  =  [1/(2**  )]^)/0Odp(p2/2m)fw(p^)exp[2ipf  M  ] 

Note:  for  f  =  0,  the  above  reduces  to  mean  kinetic  energy  of  the  system.  It  is  direct  to 
demonstrate  that 

(18)  E(x  +  f  pc-f  )  =  -(ft  2/8m)pf  f 

The  expectation  value  of  the  energy  is  the  diagonal  component  of  the  energy  density  matrix, 
and  yields: 

(19)  E(xpc)  =  [pkbT/2[[l-(X  2/6)(inp)xx] 

=  [kbT/2  +  wq]p 

The  term  wq  was  first  introduced  by  Wollard  et  al  (1990).  Equation  (19)  is  valid  only  to 
second  ordef  in  ft.  Thus  the  density  multiplying  the  term  \Jtnp)^  is  strictly  only  the 
classical  density.  The  significance  of  equation  (11)  is  that  there  is  a  change  in  energy  due  to 
the  gradients  of  carrier  density,  as  first  predicted  by  Wigner  (1932). 

THE  APPROXIMATE  NONEQUILIBRIUM  DENSITY  MATRIX 

Classical  moment  balance  equations  are  often  obtained  through  representing  the 
nonequilibrium  state  by  a  displaced  Maxwellian,  exp-[^|(p-p(j)2/2m+ V}]  where,  e.g.,  p^, 
the  density  and  a  particle  temperature,  are  to  be  determined.  The  Weyl  transformation 
indicates  that  a  displaced  Maxwellian  yields  a  density  matrix  with  the  following 
modification: 

(20)  p  (x  +  f  pc-f  y>p  (x  +  f  pcT  )exp-[2ipdC  M  ] 

Thus  the  classical  density  matrix  corresponding  to  the  displaced  Maxwellian  is  (21) 
p  (x  +  c  M )  =  Nexp-[C 2  A 2  +  p V  +  2ip^C  /*  ] 

Note  that  the  density  matrix  contains  an  imaginary  part,  as  required  for  a  finite  current. 

The  nonequilibrium  quantum-corrected  density  matrix  is  given  by  a  modification  of 
equation  (11): 
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(22) 


p  (x  +  C  PK  )  =  Nexp-[C  2  A 2  +  P V  +  2ip^C  M  ] 

{ 1-(A  2£/6)(l  +  C 2  A  2  )VXX  +  A 2  2  Vx2/12} 

Quantum  Corrected  Current  (velocity  flux^  Density  Matrix:  Since  the  density  matrix 
contains  imaginary  components  we  are  interested  in  the  expectation  values  of  the  current 
flux  density.  To  obtain  this  we  define  a  current  density  matrix  as: 

(23)  j(x  +  f  M  )  =  [l/(2*ft  )L*J0D  dp(p/m)f^ppc)exp[2ipf/ft  ] 

Note:  for  f  =  0,  the  above  reduces  to  mean  velocity  flux  of  the  system.  It  is  direct  to  show 
that: 

(24)  j(x  +  f  pH  )  =  (ft  /(2mi))p  ^ 

The  expectation  value  of  current  is  the  diagonal  component: 

(25)  j()yc)  =  (pd/m)p 

The  quantum  corrections  to  the  current  are  obtained  from  the  quantum  corrections  to  the 
density.  Note  that  the  form  of  the  current  density  along  the  diagonal  is  the  same  for 
classical  transport  as  it  is  for  quantum  corrected  transport.  Thus  of  relevance  are  the 
equations  defining  p^,  density  and  temperature.  To  establish  these  we  reconsider  the 
equation  of  motion  of  the  density  matrix  taking  successive  derivatives. 


THE  APPROXIMATE  NONEQUILIBRIUM  BALANCE  EQUATIONS  AS 
OBTAINED  FROM  THE  EQUATION  OF  MOTION  OF  THE  DENSITY  MATRIX 

Particle  Balance:  The  first  balance  equation  is  obtained  by  rewriting  the  equation  of  motion 
of  the  density  matrix  and  then  dealing  only  with  the  diagonal  component.  Thus  there  is 
reduced  content.  Using  the  definition  of  velocity  flux,  equation  (9)  becomes: 

(26)  p(x +  C  PH  )]t  +  B(x  ■ +  f  pH  )]x  -(1/i*  )[2f  Vx  +  f  3  (x +  fPH) 

+ 27  f  (p  (x  +  f  pc-f  )]f  +  [4Df 2  It, 2  ]p  (x  +  f  pc-f )  =  0 

The  first  balance  equation  is  obtained  from  the  diagonal  component  of  equation  (24): 

(27)  [p]t  +  bPd/m]x  =  0 
which  is  the  equation  of  continuity. 

Momentum  Balance:  To  determine  the  next  governing  equation  to  second  order  in  h  we 
take  derivatives  of  equation  (2)  with  respect  to  T  \  multiply  by  hi 2i,  and  obtain: 

(28)  [p(x  +  f  + 2[E(x  +  f  pt-f  )]x  +  ( 1/2)[2VX  +  f 2  V^Jp  (x  +  f  pc-f  ) 

-(1/ift  )[2f  Vx  +  f3  Vxxx/3]p(x  +  f  pc-f) 

+  2-1  p(x  +  f  pc-f  )-i[4Df  If,  ]p(x  +  C  pc-f ) 

-ifi  7  f  [p  (x  +  f  pc-f  )]f  f  +  [4D r/*3  ][p(x + f  pc-f )]  =  0 
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where  p(x  +  r  ,x-C )  =  j(x  +  f  ,x-C  )m.  The  diagonal  component  of  energy  under  finite  current 
conditions  is  required  for  the  second  balance  equation.  It  is: 

(29)  E(xpc)  =  [Pd2/2m  +  (kbT/2)(l-(X  2/6)(i  np^lp 

With  the  diagonal  components  of  energy  given  by  equation  (29),  identifying  a  relaxation 
time  r  =2/ 7,  the  second  balance  equation  is  obtained  from  the  diagonal  components  of 
equation  (29): 

(30)  (pp<j)t+(pPd2/m)x+(pkT)x  +  /7(Q/3)x+pVx+pp/r  =0 

where  we  recognize  that  [(1  np  )xx/?]x  =  -4mpQxfn  2 .  The  above  equation  differs  from  its 
classical  analog  through  the  presence  of  the  quantum  potential.  The  form  of  the  scattering 
term  in  the  above  equation  also  confirms  Frensley’s  (1990)  statement  of  the  first  part  of  the 
Fokker-Planck  scattering  as  a  frictional  term. 

Energy  Balance:  If  we  regard  the  quantum  corrected  density  matrix  as  a  distribution 
function  with  three  undetermined  parameters,  the  third  being  the  temperature,  then  a  third 
equation  is  needed  to  complete  the  system  of  equations.  This  third  equation  is  obtained  in 
a  manner  similar  that  of  equation  (28),  namely  by  taking  a  second  derivative  of  the  density 
matrix  equation  of  motion,  with  respect  to  f .  We  are  only  interested  in  the  diagonal 
components  of  this  equation,  which  with  the  energy  given  by  equation  (18),  are: 

(31)  Et  +  (l/(2m2)P3x+jVx  =  (Et)coll 

where  P3  is  the  diagonal  component  of  the  third  moment  of  the  off-diagonal  element: 

(32)  [l/(2*ft  )Uf>dp(p3  )fw(p^)exp[2ipf  jh  ]  =  (*/2i)3  Pf  f  f 
Using  equation  (21)  for  the  density  matrix,  equation  (32)  becomes: 

(33)  Et  +  {(pd/m)[E  +  (p//?)(HX  2/6](Anp)xx)]}x  +  (ppd/m)Ux  +  2E/r  +  8Dp  =0 

where  the  energy  in  equation  (33)  is  given  by  equation  (29)  (Note:  the  equilibrium  value  of 
P3  is  zero).  The  above  equation  can  be  rearranged  to  read: 

(34)  Et  +  (PdE/m)x  +  (pdpkbT/m)x  (pPd/m)[Q/3  +  V]x 

-p  (X 2  kbT/6)  [(1  np  )xx](Pd/m)x + 2E/r  -8Dp  =  0 
If  we  assume  equilibrium  values  for  D,  namely  D  =  mk^T/r,  equation  (34)  becomes: 

(35)  ^  +  (pdE/m)x + (pdpkbT/m)x + (ppd/m)[Q/3  +  V]x 

-p  (X  2  kbT/6)[(*  np  )xx](Pd/m)X  +  (2/r  )[E-kbT/2]p  =  0 

And  the  second  part  of  the  Fokker-Planck  dissipation  involves  a  relaxation  to  a  non-zero 
thermal  energy. 

The  equation  of  motion  of  the  density  matrix  involves  two  NxN  matrices  (real  and 
imaginary  parts)  where  N  represents  the  numbers  of  grid  points  along  a  linear  dimension. 
The  density  matrix  equation  is  the  equation  of  interest.  An  approximate  representation  of 
the  equation  of  motion  of  the  density  matrix  is  provided  by  the  first  three  moments  of  the 
density  matrix.  These  moments  include  the  continuity  equation,  momentum  balance,  and 
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energy  balance. 


As  in  all  moment  equation  representations,  the  form  of  the  moment  equation  depends  upon 
the  primaiy  equation,  in  this  case  equation  (9),  and  the  generic  form  of  the  distribution 
function,  in  this  case  equation  (22).  The  most  dramatic  consequence  of  the  approach,  as 
expressed  in  the  momentum  balance  equations,  is  the  appearance  of  the  quantum  potential. 
The  quantum  potential  with  its  associated  factor  of  T/3\  the  latter  arising  solely  as  a 
consequence  of  the  expansion  of  the  Liouville  potential,  permits  a  quantum  mechanical 
"correction"  description  of  classical  transport  and  allows  tunneling  to  be  incorporated  into 
the  latter.  This  result  was  also  the  initial  conclusion  of  Ancona  and  Iafrate  (1990). 

It  is  important  to  reiterate  that  one  key  result  of  the  present  study  is  that  through  the 
expansion  of  Liouville  potential  in  the  equation  of  motion  of  the  density  matrix  the 
quantum  corrections  are  the  same  as  those  obtained  from  the  equation  of  motion  of  the 
Wigner  distribution  function.  The  origin  of  the  factor  of  T/3’  is  due  solely  to  the  expansion 
of  the  Liouville  potential,  and  is  not  limited  the  high  temperature  limit  discussed  by  Ancona 
and  Iafrate  (1989). 

The  energy  balance  equation  is  also  driven  by  the  quantum  potential,  but  additionally  there 
appears  to  be  a  contribution  that  may  behave  as  a  quantum  correction  to  the  pressure.  This 
is  the  terms  p  (\ 2  k^T/6)[(X  n p  )xx\. 

In  addition  to  the  quantum  mechanical  contributions,  the  moment  equations  include  an 
incoiporation  of  Fokker-Planck  dissipation  in  the  moment  formulation,  and  twenty  of  the 
relative  contribution  of  each  of  the  two  terms.  A  more  general  treatment  is  provided  by 
Stroscio  (1986). 

The  next  question  concerns  the  significance  of  the  quantum  potential. 

SELF  CONSISTENT  NUMERICAL  SOLUTIONS 

This  section  contains  a  discussion  of  the  interpretation  of  the  distribution  of  charge  within 
the  quantum  well  and  quantum  barrier  device  and  the  values  obtained  thereof.  The 
calculations  involve  solutions  to  the  density  matrix  coupled  to  Poisson’s  equation: 

(36)  j/ax[€  (x)aV/ax]  =  -e2[p(x^)-p0(x)] 

Classical  N-^N-N^h.  Structures:  Since  most  resonant  tunnel  structures  are  designed  with 
the  heterostructures  placed  within  the  interior  of  a  low  doped  region,  the  first 
self-consistent  problem  discussed  involves  transport  through  a  classic  N  +  N’N  +  structure. 
The  structure  is  1200A  long  with  a  nominal  doping  of  10 1  3 /cm3  and  a  centrally  placed 
250A,  10 1 6  /cm3  region.  The  variation  in  background  doping  was  over  one  grid  point  or  4A. 

There  are  a  variety  of  questions  to  be  addressed  here;  among  them  is  the  issue  of  whether 
there  are  any  quantum  contributions  associated  with  an  ostensibly  classical  structure.  It  is 
relevant,  in  this  matter  to  recall  one  of  the  conclusions  of  a  paper  by  Iafrate,  Grubin  and 
Ferry  (1981),  where  it  was  argued  that  quantum  corrections  are  important  if  the  density 
sustained  an  approximate  Gaussian  distribution  with  a  width  at  half  maximum  of 
approximately  80A  of  less. 

The  equilibrium  charge  density  and  potential  distribution  as  obtained  from  the  complete 
equation  of  motion  of  the  density  matrix,  equation  (81  without  dissipation,  for  the 
N  +  N'N  +  structure  are  shown  in  figures  (la)  and  (lb),  respectively.  The  results  appear 
classical.  We  also  show  the  diagonal  component  of  energy  density  of  the  system,  figure  (1c), 
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ar  computed  from  equation  (18)  and  the  quantum  potential,  figure  (Id),  as  computed  from 
equation  (1).  For  completeness  we  also  snow  the  two  dimensional  plot  of  the  density  matrix 
in  figure  2. 

Figures  (la)  and  (lb)  display  charge  density  and  self-consistent  potential  distributions  that 
are  ostensibly  classical.  There  is  an  increase  in  potential  across  the  N‘  region,  which 
accompanies  a  decrease  in  charge  density  across  this  same  region.  From  the  point  of  view 
of  quantum  transport,  we  are  also  dealing  with  a  tunneling  problem,  particularly  with  those 
carriers  whose  energy  is  below  the  potential  barrier,  which  in  this  case  has  a  height  of 
approximately  45  mev.  We  note  that  the  mean  energy  of  the  entering  carriers  is  k^T/2  <  45 
mev. 

Quantum  transport  permits  tunneling,  thereby  increasing  the  numbers  of  particles  in  the  N~ 
region  over  the  classical  value.  Continuity  of  the  wavefunction  and  its  denvative  through 
the  barrier  region  (at  the  metallurgical  interface  the  energy  density  is  approximately 
kuTp/2,  and  decreases  into  the  N"  region;  prevents  the  density  from  approaching  it's 
classical  value,  and  instead  assumes  a  smaller  value.  This  result  which  is  a  consequence  of 
wave  function  continuity  has  been  refe  Ted  to  by  others  as  quantum  "repulsion”  (Kluksdahl 
et  al  (1989)).  The  "increased"  value  ot  density  in  the  barrier  region  ana  "decreased"  value  of 
density  in  the  classically  accessible  region  are  represented  by  equation  (16)  if  the  quantum 
is  negative  within  the  barrier  region  and  positive  within  the  region  of  the  ’turning  points’  at 
the  boundary  of  the  metallurgical  junction.  The  quantum  potential  as  evaluated  from 
equation  (1)  and  depicted  in  figure  Id  displays  these  qualitative  features.  Note  that  the 
quantum  potential  at  the  center  of  the  structure  is  approximately  25%  of  the  energy  at  the 
boundaries.  For  a  very  wide  N“  region  the  density  at  the  center  region  of  the  structure  will 
be  relatively  constant  and  Qs  0,  at  the  center. 

For  the  calculation  of  figures  1  and  2,  the  density,  potential  and  quantum  potential  are 
respectively 

(37)  ^cal(v  =  °)  =  T847  x  101 7 /cm3 

(38a)  Vcal(x  =  0)  =  44.725  mev 

(38b)  Qcal(x  =  0)  =  -3.0228  mev 

Calculations  p(x  =  0)  from  equation  (16)  with  the  potential  energy  and  quantum  potential 
given  by  equation  38  yield 

(39)  ^Eq.(38)(x  =  0)  =  1.832X101  7/cm3 

The  energy  density  at  the  center  of  the  structure  as  computed  from  the  density  matrix  is 

(40)  Eca](x  =  0)  =  2.3015  mev/cm 3 
The  energy  as  computed  from  equation  (19)  is 

(41)  EEq.(i9)(x  =  0)  “  2.286  mev/cm3 

The  hydrodynamic  moment  equations  appear  in  the  case  of  the  simple  N  +  N“N  +  to  yield 
the  same  result  as  the  density  matrix  for  the  charge  distribution. 
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Double  Barrier  Diodes:  The  situation  for  double  barrier  diodes  is  quantitatively  different 
than  that  of  the  N  +  N~N  +  structure.  Here  the  quantum  mechanics  is  not  a  correction  to 
the  classical  solution.  However,  as  in  the  N  +  N‘N  +  structure,  it  is  expected  that  outside  of 
the  double  barrier  that  continuity  of  the  wave  function  would  again  yield  a  density  below 
that  of  the  classical  value.  Thus  the  quantum  potential  would  be  positive  in  this  region. 
Within  the  barrier  quantum  mechanical  tunneling  permits  a  greater  number  of  carriers  than 
those  associated  with  tl  **  classical  distribution;  thus  the  quantum  potential  is  negative. 
Within  the  quantum  well  continuity  of  the  wavefunction  results  in  a  charge  distnbution  that 
is  below  that  obtained  classically;  and  Q  is  negative.  The  variation  of  density,  potential, 
energy  and  quantum  potential  are  displayed  in  figures  3  and  4. 

Figure  3a  displays  the  density  distribution  which  shows  a  small  residual  change  in  the 
quantum  well  and  an  insignificant  amount  of  charge  within  the  barriers.  Potential  energy  is 
displayed  in  figure  3b,  where  we  note  that  the  low  density  in  the  quantum  well  region 
contributes  to  an  elevation  of  the  potential  relative  to  the  end  points.  The  mean  energy 
density  is  shown  in  figure  3c,  and  displays  negligible  values  at  tne  center  of  the  structure. 
The  quantum  potential  is  displayed  in  figure  3d.  For  the  calculations  of  figure  3,  the  density 
potential  energy  and  quantum  potential  are  respectively 


(42)  Pcal(x  =  0)  =  6.156x  101 6/cm3 
(43a)  Vcai(x  =  0)  =  15.479  mev 

(43b)  Qcai(x  =  0)  =  92.149  mev 


Calculating  p(x  =  0)  from  equation  (16)  within  the  potential  energy  and  quantum  potential 
given  by  equation  (43)  requires  that  the  quantum  ±  Jtential  be  multiplied  by  a  constant 
other  than  1/3.  For  the  present  situation  the  constant  is  closer  to  2/3.  This  difference  is 
value  is  not  surprising  in  light  of  the  fact  that  within  the  quantum  well  the  value  of  the 
quantum  potential  is  not  longer  a  correction.  Indeed  Q  is  approximately  equal  to  the 
quasi-bound  state  energy. 

The  value  of  p(x  =  0)  from  equation  16,  with  Q/3  replaced  by  2Q/3  is: 

(44)  p(pq  (38)(x  =  0)  =  5.043  x  10 1 8 /cm3 

We  note  that  the  density  upstream  of  the  emitter  is  satisfactorily  represented  by  equation 
(16).  It  would  appear  that  our  results  are  consistent  with  the  discussion  of  Ancona  and 
Iafrate  (1989)  where  the  quantum  potential  appearing  in  the  moment  equations  should  be 

(45)  Phenomenological  quantum  potential:  aQ(x) 
where  a  is  a  position  dependent  positive  constant  less  than  unity. 

The  energy  calculated  from  equation  (18)  with  the  center  of  the  quantum  well  is 

(46)  Eca](x  =  0)  =  2.864  mev/cm3 

The  energy  calculated  from  equation  (19)  with  Q/3  replaced  by  2Q/3  is: 

(47)  E]Bq.(19)(x  =  0)  “  2.684  mev/cm3 

The  low  value  is  present,  even  though  the  quantum  potential  tends  to  push  up  the  mean 
carrier  energy,  and  arises  from  the  reduced  charge  in  the  quantum  well. 
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CONCLUSION 


This  study  represents  the  first  comprehensive  effort  to  assess  through  (1)  numerical 
simulations  of  the  equation  of  motion  of  the  density  matrix,  and  (2)  approximate  analytical 
procedures  involving  an  expansion  of  the  Liouville  potential,  the  multiplicative  constants 
associated  with  the  quantum  potential,  as  used  in  quantum  hydrodynamic  transport.  It  was 
found  that  when  quantum  effects  are  ’corrections’  to  classical  calculations,  as  in  N  +  N’N  + 
structures  the  factor  T/3’  is  a  satisfactory  multiplicative  constant.  When  the  quantum 
potential  is  no  longer  a  correction,  as  in  the  case  of  double  barrier  diodes,  the  multiplicative 
constant  is  no  longer  ’constant’  but  is  position  dependent  and  less  than  unity.  Additional 
work  is  required  to  narrow  the  range  of  variation  of  this  constant  before  effectively 
incorporating  the  quantum  potential  in  detailed  quantum  hydrodynamic  simulations  of 
ultrasmall  devices. 
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Fig.  lb.  Self  consistent  potential  associates  with  the 

density  distribution  of  Figure  la.  The  N+N"  and 
N"N+  regions  are  indentified  by  the  dashed  lines 
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lc.  Energy  density  as  obtained  from  Equation  (18) . 
Density  at  boundary  region  is  entirely  thermal. 
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Fig.  2.  Two  dimensional  representation  of  the  density  matrix,  p(x,xf),  for 
the  NIN  structure  whose  diagonal  components  are  shown  in  Figure  1. 
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Fig.  3b.  Self  consistent  potential  plus  heterostructure  potential 
associated  with  the  charge  distribution  of  Figure  3a.  The 
N+N~  and  N“N+  regions  are  identified  by  the  dashed  lines. 
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Quantum  potential  for  the  charge  distribution 
of  Figure  3a  as  obtained  from  Equation  (1) . 
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Fig.  4.  Two  dimensional  representation  of  the  density 

matrix,  p(xfx*)f  for  the  double  barrier  structure 
whose  diagonal  components  are  shown  in  Figure  3. 


